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' fj , Abstract- An evolution operator £„ with n arbitrary, typical of several models, is analyzed. When 

tZJ I n = 1 the operator characterizes the Standard Linear Solid of viscoelasticity, whose properties are 

already extablished in previous papers. The fundamental solution £„ of £.„ is explictly obtained and it's 
estimated in terms of the fundamental solution £i of £i. So, whatever n may be, asymptotic properties 
and maximum theorems are achieved. These results are applied to the Rouse model and reptation model, 
jrt ' which describe different aspects of polymer chains. 
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1.INTRODUCTION 

^ I Let Cn be the (2 + n) order operator: 

■^ ' 

p 

g : Cn = #^(öit - cld^^) + a„_i d^;'-^\du - cLiö.x) + •■ (1.1) 

(N 



•••OO (dtt - CqÖ^x) 



^ I where a^ {k = 0..n — 1) are positive constants. 



According to the value of n, (1.1) describes several physical phenomena. As example, 
when n = 1, Cn can be found in dynamic of relaxing gases, in magnetohydrodynamics, 
in hereditary electromagnetism (see [1] and references therein) and in isotropic viscoelas- 
ticity where, (1.1) models the evolution of the Standard Linear Solid (S.L.S.) (see, f.i. 
[2, 3]). 

In all these models c^ 's represent the characterized speeds depending on the materials 
properties of the medium and in many physical problems it results Cq < Ci-. < c^_i < c^ 
as it's typical of wave hierarchies. [4]. 

When n = 1, the operator (1.1) is stricty-hyperbolic and it has been widely anal- 
ized in [1]. It's fundamental solution £i has been explicitly determined and singular 
perturbation problems, together with asymptotic properties, have been estimated. 

Aim of the paper is to draw generalizations of the wide analysis related to S.L.S. 
to the case of (1.1) with n arbitrary. For this, a conditioned equivalence between (1.1) 
and an integro - differential operator Ai related to an appropriate memory function 
gn{t), is considered. Owing to the hypotheses of fading memory, every function gn{t) 
can be approximated by Dirichlet polynomials [5, 6] with appropriate restrictions on the 



coifïicients of this expansion. These limitations need that the differential operator £„ is 
typical of wave hierarchies [7, 8]. 

By this equivalence, whatever n may be, the fundamental solution Sn of (1.1) is 
expUctly achieved and it is estimated in terms of £i. So the niaxiniuni properties and 
asymptotic estimates established for the Standard Linear Sohd can be apphed to operator 
Cn defined in (1.1). Moreover, boundary layer problems, typical of dissipative media, 
can be rigorously estimated. [9, 10]. 

These results are apphed to the Rouse model and the reptation model which describe 
different aspects of polymer chains and have met with reasonable success.[ll]-[13]. 



2. DIFFERENTIAL CONSTITUTIVE EQUATION 



Let B a hnear, isotropic, homogeneous system and let u{x, t) the displacement field 
from an undeformed reference configuration Bq. If u = u{x,t) i and po denotes the mass 
density in Bq, the equations of one-dimensional motions of B are 

PqUu = (7 + f, e = u, (2.1) 

where ƒ = ƒ! is the known body force, while a and e are the only non vanishing 
components of the stress and the strain tensors. 

When the viscoelastic behavior of B is of rate-type, the well known stress-strain 
constitutive relation is 



fc=0 fc=0 

with ak,ak constant {an,an 7^ 0). 

Then, by (2.1), (2.2), the displacement field u{x,t) is solution of the higher order 
equation like: 



^nV = ^ ak dt {vtt - CkV^x) = F (2.3) 

A;=0 

where: 



Cfc = ak/poük, F = (l/po) Yl "fe ^tf- (2-4) 

fc=0 

The constitutive relation (2.2) includes various classical mechanical models, as Maxwell 
and Keivin - Voigt models [14]. Moreover, when n = 1 and 

O < Co < Cl, 7] = ai/ao > O (2.5) 



one has the case of the Standard- Linear Solid (S.L.S.) which is niodelled by the strictly- 
hyperbohc third order equation: 

Civ = 7] dt{vtt - ciVxx) + Vu - Co v^x = {l/ao)F. (2.6) 

The fundamental solution £i of this operator has been obtained in [1], also when 
X £ R^ OT X £ R^. Further, numerous basic properties of £i have been rigorously 
estimated and the wave behavior of S.L.S. is now acquired. 

When n > 1 and all the Ck s are positive, then waves of different orders appear and 
their roles must be clarified in order to see how each set is modified by the presence of 
the other. Obviously, wave or dispersive behavior depend on the requirements of the 
coefïicients a^ and Ck due to physical properties of the system B. 

For this, we will analyze the restrictions imposed on the constants a^ and c^ by usual 
hypotheses of fading memory for B. 



3.INTEGRAL CONSTITUTIVE EQUATION 



When the strain amplitudes are not too large, the behavior of most viscoelastic media 
is fairly well modelled by linear hereditary equations like: 

e{t) = J(0)[cj(t) + / j{t- T)aiT)dT], (3.1) 

where J{t) denotes the creep-compliance and the integral term needs the knowledge of 
the past history of the stress. 

Usually, according to fading memory hypotheses, J{t) is a positive fast decreasing 
function. For instance, several real materials as polymers, rubbers, bitumines, have 
satisfactory representations by means of chains of S.L.S. elements in series or parallel 
[3, 11]. In the series case, the creep function is 

n r> 

J^{t) = U0)[l + J2-±il-e-^''% (3.2) 

where n is the number of elements in the chain, tj. = /3^ are the characteristic times 
and Jn(0) denotes the elastic compliances. 
Then, if one puts: 

c' = [poJniO)r\ F, = c2[J„(0)/ + / Jnit - T)(7,(T)dT], (3.3) 

by (2.1), (3.1), (3.2) one deduces: 



Mu = c^Uxx - utt - g{t-T)urTdT = -F^{x,t), (3.4) 



with 



n 



9 = 9n{t) = Y.Bke-^''' =Jn{t)/Jn{0). (3.5) 



fc=l 



In this memory function, n is quite arbitrary and constants Bk and frequencies f3k are 
such that: 



O < /3i < /?2... < /3n; Bk>0 V/c = 1, 2....n. (3.6) 

These hypotheses assure that: 

g{t)>0, g<0, g>0, Vt > O (3.7) 

according to the convexity assumption considered by Dafermos [15]. 

We observe that the representation (3.5) of the memory function g is not restrictive 
because well-known Muntz and Schawart'z theorems [5, 6] imply that whatever C^{R^) 
function can be uniformly represented by means of Dirichlet polynomials. Moreover, as 
n is arbitrary, the constants -Bfc,/3fc can be determined in order to fit the experimental 
curves for g{t) to any prefixed degree of approximation [11]. 

By (3.2), (3.6) one has 

n Tj 

J„(oo) = J„(0)[1 + V -i] > J„(0). (3.8) 

k=i « 

4. FADING MEMORY AND WAVE HIERARCHIES 

The initial data related to (2.3) and (3.4) let be null and let 



p{s) = E /"fc ^' ' Qi') = E ^-fc ^' (4-1) 

k=0 fc=0 

with 



/"fc = ak/ttn , Afc = akCk/ünCn {k = o, ..n) (4.2) 



so that iJ-n = Ki = 1- 
Further, let 



n o 

the Laplace transform of the memory function (3.5). 

Then, if one apphes the Laplace transformation to (2.3) and (3.4), it results: 

c„ (;>(s) a„c„ Q(s) 

üxx-^ [l + Gis)]ü= -J (4.5) 

where (') denote the L -transform of ( ). 
By comparing (4.4), (4.5) one deduces 

and the pohnomial identity imphes c„ = c^ and 



Ao = /3i/32../3„ 

An-2 = /3i/32 + /3i/33 + ../3„-i/3„ 
A„_i = /3i + .... + /3„ 

So, owing to (3.6), all the A^'s are positive. Further, as for n^, one has: 
^0 = Ao + 5i(/32../3„) + ...B„(/3i../3„_i) 



(4.7) 



(4.8) 



^n~2 = A„_2 + Bi{f32 + .. + /3„) + ...Bn{f3l + .. + /3„_i) 
^n-1 = A„_i + Bi + ....+ Bn 

and (3.6), (4.8) imply too: O < Afe < /ifc (/e = O, ..n — 1). As consequence: 



O < Cfc < c„ = c^ (A: = O, ...n - 1) (4.9) 
At last, by (4.7), (4.8), it follows: ^ < ^ < ^^ and so 

O < Co < Cl.... < c„. (4.10) 

So, the foUowing property holds: 

Property 4.1. Hypotheses of fading memory (3.5) (3.6) ivnply that the differential 

operator (2.3) is typical of wave hierarchies. I 



Vice versa, the inverse transformation of (4. 7), (4. 8) requires carefulness. 
When the difFerential equation (2.3) is prefixed, in order to obtain the dual hereditary 
equation (3.4) with a memory function g{t) satisfying (3.5), (3.6), appropriate restric- 
tiones on the constants ak,Ck must be imposed. 

At first, (4.3), (4.6) imply 

Pis) ^ Bk , , 

-f4 = Bo + Y] — 4.11 

Qis) t^,s + jk 

where all the roots s = — 7^ di Q{s) are real and simple, with 7^ > 0. Moreover the 
conditiones B^ > O, which are sufficiënt to verify (3.7), involve further hmitations. 
Example 4.1 - When n = 1, one has: cP = ci, B^ = 1, and 

/3i = ^>0 ,i?i = ^(l-^)>0 (4.12) 

aici ai Cl 

which represent the known restrictions typical of S.L.S. I 

Example 4.2 - When n = 2, one has: cP = C2, Bq = 1, and /3i, /32 are real iff: 

u;2 = ^aicif - 4(aoCo)(a2C2) > O (4.13) 

Then, it results: 

/3i = (aici - w), 132 = (aici + u) (4.14) 

2a2C2 2a2C2 

so that O < /3i < /32 . Further 

Bi = ^ [ao(c2-co)-aiA(c2-ci)]. (i = l,2) (4.15) 

ÜJ 



Thus, it is Bi > O, 02 > O iff 



/3i < ^ ^^^ < /32. (4.16) 

öl C2 - Cl 



Therefore, the fourth-order operator 



a2{utt - C2Uxx)tt + ai{utt - ciUxx)t + ao{utt - cqu^x) (4-17) 

can be analyzed by (3.4), (3.5) , (3.6) when the constants ak,Ck satisfy (4.13) and (4.16). 

■ 



5. ESTIMATES FOR THE HEREDITARY MODEL 

Let Bn the viscoelastic model characterized by the memory function gn in (3.5); the 
case n = 1 corresponds to the L.S.L. Bi. 

In [16, 17], the fundamental solution En of the operator Ai in (3.4) has been expUcitly 
determined, whatever n may be. If rj is the step - function and /q is the modified Bessel 
function of first kind, it results: 

En = En{l5i..pn, ^l-^n) = 7^r?(t - r){Ai + A2) (5.1) 

2c 

with 

Al = e-^«*/2^o(f Vt^^72) (5.2) 

A2 = - r d9 f e"ö"^//(z, t - w)du, (5.3) 

TT Jo Jr 

and go = g{0),r = \x\/c, 2z = w — cosO {nP' — r^Y ■ Futher, if 

,/.fc(z,t) = e-^'^'^Bkhz/t h{2^Bkhz t), (5.4) 

one has: 



H{z,t) = Y^ (l)k+ Yl (Pki * 4>k2 + - (5.5) 

ki,K2 

where sums are computed according to the simple combination of the indices kik2--kn 
and * denotes the convolution with respect to t. 

Moreover the fundamental solution En related to Bn and defined in (5.1)-(5.5), can be 
rigorously estimated in terms of the fundamental solution Ei related to an appropriate 
S.L.S. BI defined by 

n D 

5i=6e-^i* yjHh b = ^lY.-^. (5.6) 

l Pk 

In fact, if r is the open forward characteristic cone {{t,x) : t > O |x| < et}, and 
Xn = 11^=2 ("5^)^5 then the following theorem holds: 

Theorem 5.1 - If the memory function is given by (3.5) (3.6), then the fundamental 
solution En of M. is a never negative C°°{T) function and it satisfies the estiniate: 



0<Enif3l../3n,Bi..Bn)<Xn Ei{^i,b), (5.7) 



everywhere in T and whatever n niay he. I 

Remark 5.1 - The model B\ defined by (5.6) is physically meaningful. 
In fact the memory function gi is related just to the obhviator because ri = /3j~ is the 
longest characteristic time. 

Furthermore ö„ and B\ verify the same hypotheses of fading memory and by (3.5) 
(5.6), it results: 

/ gn{t)dt = / gi{t)dt =Y.ir- (5-8) 

Jo Jo ]^ Pk 

Moreover, the integral (5.8) affects the asymptotic analysis of hereditary equation [7, 18]. 

■ 

Remark 5.2. - By known properties of asymptotic behaviour of convolutions, the 
constitutive relation (3.1) imphes: 

hm e(x,t) = J„(0) [1+ / gn(t)dt] lim a{x,t) (5.9) 

provided that linii^^oo o" exists. Then Bn and the model Bf exhibit the same asymptotic 
behaviour (5.9) 

Further, by (3.5), (3.8) it results: 

and so, (5.10) implies J„(oo)(t(oo) = e(oo). Consequently, when t is large, the behaviour 
of Bn is typical of an elastic material with modulus J„(cx)). I 



6. ESTIMATES RELATED TO WAVE HIERARCHIES 



When the operator £„ is reduced to the hereditary operator A4 of (3.4), then esti- 
mates of Theorem 5.1 can be applied to wave hierarchies. Obviously, the equivalence is 
conditioned by inverse transformation of (4.7)- (4.8) together with (3.6) (see n.4). 

Let Cl the operator (2.6) related to the S.L.S. Bf characterized by 

Pi+b Pi Pi + b 

Now, let Cn the differential operator given by (2.3) whatever n may be, and let Vn 
a prefixed boundary value problem related to £„• The meaningful aspects of qualitative 
analysis of the solution of Vn can be obtained by means of Theorem 5.1 and by the 
known properties of Cl [1]. 



So maximum properties, asymptotic behaviour, boundary layer estimates, etc. for 
the solution of P„ are deduced by analogous properties related to £ J . 

Moreover, owing to the equivalence between £„ and A^, it is possible to have exphctly 
the fundamental solution of Cn for all n. In fact it sufïices to apply the explicit formula 
(5.1)-(5.3). 

As an example the case of polymeric materials can be considered. 

Example 6.1 - Polymeric materials are very flexible and are easily formed into 
fibres, thin films, additives for oils, etc. So theirs applications to concrete problems are 
numerous.[12, 19, 20]. According to theories of linear viscoelasticity, two models, that 
descibe different aspects of polymer chains, have met a reasonable success: the Rouse 
model and the reptation model [13]. 

In both cases the memory function g{t) assumes a form likc (3.2) (3.3). In fact in the 
reptation m,odel, the stress relaxation function is: 



where k is a constant depending on the polymer physics and the value of the "reptation" 
time Til can be fixed according to elasticity experiments [11]. 

When the viscoelastic behaviour is represented by the Rouse model, memory function 
g{t) is given by : 



ö(t) = A;i^e'"'n, (6.3) 

h=l 

where the relaxation time ri can be calculated by means of experimental results.[13]. 

So, if one considers the first two steps in the reptation model, it results: Bi = k, B2 = 
Bi/9, /3i = 1/Td, P2 = 9/3i. Consequently the operator (4.17) is characterized by 
constants: 



Cn = c^ ^ Cl = C^ 9 C9 = c^ 

'-o <- 81+82fcrd 1 ^ 9+fcrd ^2 t, 

(6.4) 



Analogously, in the Rouse model, beeing Bi = B2 = ki, /3i = 2/ri, /32 = 4/3i, one has: 



co = C' 8+^777 Cl =0^5^ C2=c2 



^2 _,„ -2 



(6.5) 



ao = l + ^ri ai = ^(2A;i + i^) 03 = ^. 



The wave hierarchies defined by (6.4) or (6.5) are governed by the operator C^ of the 
Standard Linear Solid defined, respectively, by: 



^0 — '- 81+82fcrd '-1 — ^ "O — J- ^ 81 '*''d '/ ~ 81+82fcTd ' 

(6.6) 

.co = c2g^^ ci=c2 ao = l + |A;iri r/=g4^. 

■ 
These results have been confirmed also in [21] for entangled polymers with chain 
strecth. 
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